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Abstract 



Results from the Vr appearance search in a neutrino beam using the full NOMAD 
data sample are reported. A new analysis unifies all the hadronic r decays, signifi- 
cantly improving the overall sensitivity of the experiment to oscillations. The "blind 
analysis" of all topologies yields no evidence for an oscillation signal. In the two- 
family oscillation scenario, this sets a 90% C.L. allowed region in the sin^ 20^,-— Am^ 
plane which includes sin^ 20^^^ < 3.3 x 10~^ at large Am^ and Am^ < 0.7 eV^/c^ 
at s\v?26f^r = 1. The corresponding contour in the Ue t^t oscillation hypoth- 
esis results in sin^ 2^eT < 1-5 x 10~^ at large Am^ and Am^ < 5.9 eV^/c'' at 
sin^20er = 1- We also derive limits on effective couplings of the r lepton to u^ov 

Key words: appearance, neiitrino mass, neutrino oseillations 
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1 Introduction 



The NOMAD experiment was designed in 1991 to search for i/j- appearance 
from neutrino oscillations in the CERN wide-band neutrino beam produced by 

the 450 GeV proton synchrotron. Given the neutrino energy spectrum and the 
distance between the neutrino source and the detector (~ 600 m on average), 
the experiment is sensitive to differences of mass squares Am^ > 1 eV^. 

The experiment was motivated by theoretical arguments suggesting that the 
may have a mass of 1 eV, or larger, and therefore could be the main con- 
stituent of the dark matter in the universe. This suggestion was based on two 
assumptions: 

• the interpretation of the solar neutrino deficit [1] in terms of i^e I'lj, oscil- 
lations amplified by matter effects [2], giving Am^ ^ 10~^ eV^; 

• the so-called "see-saw" model [3] which predicts that neutrino masses are 
proportional to the square of the mass of the charged lepton, or of the charge 
2/3 quark of the same family. 

From these two assumptions one expects a z/^mass of ~ 3 x 10~^ eV and a 
Uj- mass of ~ 1 eV, or higher. Furthermore, in analogy with quark mixing, 
neutrino mixing angles were expected to be small. 

It is within this theoretical scenario that the NOMAD experiment [4] started a 
direct search for i/^ Vt oscillations together with another experiment (CHO- 
RUS) which used the same neutrino beam but searched for appearance with 
a complementary technique [5] . 

The detection of an oscillation signal in NOMAD relies on the identification 
of Vt- charged- current (CC) interactions using kinematic criteria. The spatial 
resolution of the detector does not resolve the r decay vertex from the v-r CC 
interaction. The identification of Vt CC events is thus achieved by exploit- 
ing all the kinematic constraints which can be constructed from a precise 
measurement of all visible final-state particles. This requires a detector with 
good energy and momentum resolution and sophisticated analysis schemes. 
The NOMAD experiment is the first neutrino oscillation search to use this 
technique [6]. 

From 1995 to 1998 the experiment has collected 1 040 000 events with an iden- 
tified muon, corresponding to about 1 350 000 CC interactions, given the 
combined trigger, vertex identification, and muon detection efficiency of 77%. 

A recent paper [7] described a search for and Vf. oscillations in 

the full NOMAD data sample. The analysis was based on both deep inelastic 
(DIS) interactions and low-multiphcity (LM) events for all the accessible r 
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decay channels. No evidence for oscillations was observed. 

In this paper we report an improved appearance search in the hadronic r 
decay channels. We then combine this search with the analyses of the VT^Ve 
DIS channel and of the LM topologies described in Ref. [7] and we present 
the overall NOMAD results. The new analysis unifies all the hadronic r de- 
cay topologies (z/t-/i, and z/^3/i) in a single selection scheme. A refined 
implementation of kinematics, exploiting all the available degrees of freedom, 
together with new algorithms for the rejection of CC backgrounds leads to an 
improvement of the combined NOMAD sensitivity to oscillations by almost a 
factor of two with respect to Ref. [7]. The systematic uncertainties are also 
substantially reduced. 



2 NOMAD detector 



The NOMAD detector (Figure 1) is described in detail in Ref. [4]. Inside a 0.4 

T magnetic field there is an active target consisting of drift chambers (DC) [8] 
with a fiducial mass of about 2.7 tons and a low average density (0.1 g/cm^). 
The target, 405 cm long and corresponding to about one radiation length, is 
followed by a transition radiation detector (TRD) [9] for electron identification, 
a preshower detector (PRS), and a high resolution lead-glass electromagnetic 
calorimeter (ECAL) [10]. A hadron calorimeter (HCAL) and two stations of 
drift chambers for muon detection are located just after the magnet coil. The 
detector is designed to identify leptons and to measure muons, pions, electrons 
and photons with comparable resolutions. Momenta are measured in the DC 
with a resolution: 

CTp 0.05 0.008 X p[GeV/c] 



where L is the track length and p is the momentum. The energy of electro- 
magnetic showers, E, is measured in the ECAL with a resolution: 

E . ^E[GeV] 



The neutrino interaction trigger [11] consists of a coincidence between signals 
from two planes of counters located after the active target, in the absence of 
a signal from a large area veto system in front of the NOMAD detector. 
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Fig. 1. Side view of the NOMAD detector. 



3 Neutrino beam and event samples 



From recent beam computations [12], tlie relative composition of CC events 
in NOMAD is estimated to be z/^ CC : z/^ CC : CC : z/g CC= 1.00 : 0.0227 : 
0.0154: 0.0016, with average neutrino energies of 45.4, 40.8, 57.5, and 51.5 
GeV, respectively. The prompt v^- component is negligible [13]. Neutrinos are 
produced at an average distance of 625 m from the detector. In addition to 
the z/^ CC events, the data contain about 34 000 z/^ CC, 21 000 z/^ CC, 2300 
z/g CC interactions, and 485 000 neutral current (NC) interactions. 

Neutrino interactions are simulated using a modified versions of LEPTO 6.1 
[14] and JETSET 7.4 [15] with and cutoff parameters removed, and 
with r mass and polarization effects included. We use the nucleon Fermi mo- 
tion distribution of Ref. [16], truncated at 1 GeV/c. To define the parton 
content of the nucleon for the cross-section calculation we use the GRV-HO 
parametrization [17] of the parton density functions, available in PDFLIB [18]. 
A full detector simulation based on GEANT [19] is performed. Further cor- 
rections to these samples are applied using the data themselves, as described 
in Section 4.2. The size of the simulated samples exceeds the data sample by 
a factor of about 3 for z/^ CC interactions, 10 for NC and z/^ CC interactions 
and 100 for z/g and Vf. CC interactions. In addition, more than 500 000 Vr CC 
events have been generated for each r~ decay channel. 
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Fig. 2. Signal and background topologies in NOMAD: a) NC background; b) CC 
signal with subsequent r decay; c) f^(fe) CC background. The square indicates the 
reconstructed "primary" vertex for u^- CC interactions. The effect of the Ty selection 
on i'i^{i>e) CC topologies is discussed in Sections 5.3 and 5.4. 

4 Analysis principles 

Prom the kinematical point of view, v-,- CC events in NOMAD are fully char- 
acterized by the (undetected) decay of the primary r. The presence of visible 
secondary r decay products, ry , distinguishes them from NC interactions, 
whereas the emission of one(two) ncutrino(s) in hadronic(leptonic) r decays 
provides discrimination against u^^u^ ) CC interactions (Figure 2). Conse- 
quently, in Vt CC events the transverse component of the total visible mo- 
mentum and the variables describing the visible decay products have different 
absolute values and different correlations with the remaining hadronic system, 
H ^ than in [v^. ) CC and NC interactions. The optimal separation between 
signal and background is achieved when all the degrees of freedom of the event 
kinematics (and their correlations) are exploited. 

A rejection power against backgrounds of (9(10^) is required from the kine- 
matic analysis in order to match the data sample size (Section 3). In addition, 
the potential z/^ signal allowed by limits from previous experiments [20,21] 
is at least by a factor of 0.0025 times smaller than the main v^CC com- 
ponent. Therefore, the Vj. appearance search in NOMAD is a search for rare 
events within large background samples. This imposes severe constraints on 
the analysis techniques. In order to obtain reliable background estimates we 
have developed methods to correct Monte Carlo (MC) predictions with exper- 
imental data and we have defined appropriate control samples to check our 
predictions. 

This paper describes a new search for Vj. CC interactions in the hadronic r 
decay channels h~{mr'^)i'r and h~h'^h'' (n7r°)z/,- , for a total branching ratio 
of 64.7% [25]. The analysis focuses on DIS events, defined by a cut on the total 
hadronic momentum recoiling against the visible r decay product(s), > 1.5 
GeV/c. 

Neutrino interactions in the active target are selected by requiring the presence 
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of at least one charged track in addition to the potential r decay products, 
originating from a common vertex in the detector fiducial volume. Quality cuts 
are then applied to ensure that the selected events are properly reconstructed. 
This is obtained by imposing constraints on the statistical significance of the 
main kinematic variables (see Appendix A and Section 5.5): cr(p^)/p^ (trans- 
verse plane), a{QLep )/Qlcp , cr{p^ )/p^ and (T^^^{p)lp^ ([p"^ , p^] plane). In 
addition, loose requirements based on approximate charge balance at the pri- 
mary vertex are also applied. Overall, these quality cuts typically remove 10 
to 15% of the events. 

The separation of Vj. CC interactions from backgrounds is described in detail 
in Section 5 and can be summarized in the following steps (Figure A.l): 

• the visible r decay products ry are identified on the basis of their topology 
(Section 5.1); 

• constraints are applied to the structure of the hadronic system H recoiling 
against ry (Section 5.2); 

• specific algorithms are used to identify and veto primary leptons originating 
from Vfjt {Vf. ) CC interactions (Section 5.3); 

• constraints are imposed on the internal structure of the selected Ty candi- 
date (Section 5.4); 

• the final background rejection is achieved by exploiting all the available 
information from the global event topology (Section 5.5). 

Section 6 is devoted to the background estimate, with a description of the data 
corrections (Section 6.1) and of the control samples used to check background 
predictions (Section 6.2). Systematic errors are discussed in Section 7. The 
final results from the analysis of the hadronic DIS channels are then combined 
in Section 8 with the remaining topologies from Ref . [7] . 

In the following we describe the main analysis principles, which have been 
extensively discussed in Ref. [7] [26]. 

4..I Statistical analysis of data 

The kinematic variables used in the present analysis are defined in Appen- 
dix A. In order to exploit their correlations, these variables are combined 
into likelihood functions. These functions are partial two-, three- or four- 
dimensional (2D, 3D, 4D) probability density functions {pdf), denoted in the 
following by square brackets (e.g. [a, 6] denotes the 2D correlation between 
the variables a and h) . The final likelihood function, £, is then obtained from 
a combination of partial pdf 's which includes correlations (denoted again by 
square brackets) among them. This analysis technique provides the best sen- 
sitivity to oscillations. As is common practice, the logarithm of the final like- 
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lihood ratio between test hypotheses, In A, is used. 

A fit to In A with separate signal and background components can increase the 
statistical power of the search. However, since few events are expected in the 
signal region, we combine regions of similar signal to background ratio into 
suitably chosen bins of variable size (Section 5.5). These bins are then treated 
as statistically independent (Section 8.3). 

4-2 Data simulator 

In order to rehably compute signal and background efficiencies, the MC re- 
sults are corrected using the data themselves (Section 6.1). We perform this 
correction by using a sample of z/^ CC events from the data, removing the 
identified muon, and replacing it with a MC-gencratcd Icpton i of appropriate 
momentum vector, where i can be a u, an e~ or a t~ followed by its decay. 
In addition, as explained in Section 6.1, events with identified muons or elec- 
trons are directly used to correct CC background predictions in the hadronic 
T decays. All these samples are referred to as the Data Simulator (DS). 

The same procedure is then applied to reconstructed MC CC events and 
these event samples are referred to as the Monte Carlo Simulator (MCS). In 
order to reduce systematic uncertainties, signal and background efficiencies e 
are then obtained from the relation: 

e = CMC X Cos/ CMCS 

which implies that efficiencies for lepton reconstruction are obtained from 
the MC (eMc)i while the effect of the hadronic jet differences between data 
and simulations is taken into account through the ratio eDs/^MCs- It has been 
checked that efficiencies obtained by this method are indeed stable with respect 
to variations of the MC input models (nuclear effects, fragmentation, Fermi 
motion, detector resolution functions) within the quoted uncertainty. 

The errors on background predictions given in all the following tables reflect 
the statistical uncertainties from MC, MCS and DS. Systematic uncertainties 
are discussed in Section 7. 

4-3 Avoiding biases 

A procedure referred to as "blind analysis" is used. According to this proce- 
dure, data events inside the signal region (the "blind box" ) cannot be analyzed 
until all the selection criteria are defined and the robustness of the background 
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predictions is demonstrated with appropriate control samples (Section 6.2). 
The selection criteria are chosen by optimizing the sensitivity to oscillations 
(Section 5.5). This is defined as the average upper limit on the oscillation 
probability that would be obtained, in the absence of a signal, by an ensemble 
of experiments with the same expected background [27]. 

The sensitivity to oscillations also provides the final criterion for any choice 
between different analyses of the same r decay topology or for the replacement 
of a previous analysis by a newer one. The choice is made before looking at 
data events falling in the signal region. 

The MC samples used to define the selection criteria are different from those 
used to evaluate the background and signal efficiency, which are therefore fully 
unbiased. 
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5 Selection scheme 

The new selection of hadronic r decays unifies four specific decay topologies 
in a single general scheme: 

O7: T Ut-tt decays. The ry candidate is built from a single charged track. 
I7: r — >• UtP ~^ i^TTTTT^ decays where the 7r° is reconstructed as a single ECAL 

cluster either because the two photons overlap in ECAL or because one 

of them is not reconstructed in ECAL. The ry candidate is built from a 

charged track and from a single ECAL neutral cluster. 
27: r VrP ~^ VrT^ii^ decays where the 7r° is reconstructed from two separate 

ECAL clusters. The t\ candidate is built from a charged track and from 

two ECAL neutral clusters. 
3/i: T — > i^T-oi — > VrT^p^ — > Vrnnn dccays. The ry candidate is built from three 

charged tracks. 

Each topology is independently analyzed with the same selection criteria. 
Events selected by more than one topology define a further sample. Eventu- 
ally, we combine the different topologies into a single search. This statistical 
treatment provides the best overall sensitivity to oscillations. Photon conver- 
sions in the DC volume are not used for tt" reconstruction in this analysis. 
However, events where only one of the 7's from 7r° decay converts into an 
e^e~ pair are included in the I7 topology. The signal efficiencies quoted in 
the following refer to the i/^/i(n7r°) (Br 49.5%) or i/.^3/i(n7r°) (Br 15.2%) [25] 
inclusive DIS samples. 

The selection scheme is intended to exploit, at each step, all the available 
topological information (degrees of freedom) through the use of appropriate 
probability density functions based on correlations among kinematic variables. 
As explained in Appendix A, an event in NOMAD can be fully described by 
five degrees of freedom (Figure A.l): three in the transverse plane (x, y) and 
two along the beam direction {z). 



5.1 Identification of T decay product(s) 

In hadronic r decays the selection of the visible decay product (s) relies on 
topological constraints. This implies that, for a given event, more than one 
choice is possible and therefore a discriminating criterion is needed. Since the 
total visible momentum of the event (i?vis and p^) is fixed in the problem, only 
three degrees of freedom related to the general event structure are available. In 
addition, when applicable, the internal structures of ry and H provide further 
constraints. 
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A single likelihood function, , is used to select among all the combinations 
the most likely Ty originating from a r decay: 



£ = [ (£ ),RQT.,yBj,OTvH] 

The 2D correlations between the three variables Rqj,, ysj and 9ryH (Ap- 
pendix A) are shown in Figure 3 (left plots). The function £™ describes the 
internal vy structure for the I7, 2j and 3h topologies (obviously, it cannot be 
defined for the O7 topology): 



'IN 



[ Mp, 9^-^o, E^o/E^is ] I7 
[[ M,o , e^^, E^'^-IE,,, ] , Mp, ^,-,0 , E,./E,,, ] 27 

[[ MpO , ^^+7r- , Et^+ ] , , 6'7r-7r- , E pO ] 3/l 

where, for the 3/i topology, the average of the two indistinguishable tt+tt" 
combinations {M pO^O^^+T^- ^ E po) is used. This function is based on the decay 
kinematics of the p (I7 and 27) and ai (3h) resonances in the laboratory frame. 
In the definition of C^^ the invariant masses are considered as free parameters. 

The isolation variable Rq^ is sensitive to the internal structure of the hadronic 
system and thus incorporates additional information with respect to the de- 
grees of freedom related to the global event topology. The use of this variable 
for the Tv identification results, overall, in four degrees of freedom in addition 
to the Tv internal structure. However, only three of them are physically rele- 
vant for the Tv identification and are therefore included in . This procedure 
(and the choice of the variables entering ) also minimizes potential biases 
induced on background events, in which a fake (not coming from r decay) 
Tv is indeed constructed by the tv identification algorithm itself. 

A likelihood ratio. In , is built as the ratio of the function for correct 
and wrong combinations of r decay products in Vr CC events. In each event, 
the system of particles with maximum In among all possible combinations 
is tagged as tv . No cut is applied on In A^ , which is constructed from signal 
events only and, therefore, is not the optimal discriminant against backgrounds 
(Section 5.4). After the ry identification, a minimum 7 energy of 200 MeV 
(100 MeV) is required for the 17(27) topology in order to remove fake photons. 
The algorithm correctly identifies the visible t decay product(s) in 91%, 80%, 
80% and 73% of all the O7, I7, 27 and Zh v-j. CC events respectively. These 
values set an approximate upper limit to the final selection efficiency, since 
events with a mis-identified tv do not have the kinematics expected from CC 
events (Section 5.5). 
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Fig. 3. Left plots: correlations between kinematic variables, not related to the internal ry structure, used to construct In for wrong 
(a,c,e) and correct (b,d,f) combinations of r decay products in Ut-tt decays (Section 5.1). Right plots: correlations between liinematic 
variables used to construct In for wrong tracks (a,c,e) and for the true leading muon (b,d,f ) in CC events passing the first level 
veto (Section 5.3). 



5.2 Structure of the H system 



For background events the ry candidate is mostly selected inside the hadronic 
jet. In NC interactions the whole visible event is indeed a hadronic jet (Fi- 
gure 2a). On the other hand, in CC interactions where the leading lepton is 
not identified and the rvparticle(s) are embedded in the jet, the remaining 
hadronic system H has an incorrect topology (Figure 4b). As a consequence, 
a constraint on the structure of the hadronic system significantly increases the 
background rejection. 

In order to effectively reject both NC and CC topologies, the jet structure 
must include information from ry . This is achieved by defining the variable: 



{Q't)h 



which takes into account the transverse size of the hadronic system, {Q^)h, 
and the opening of the minimum invariant cone between Ty and any other 
charged track, Ar^^y/j.. The requirement Sh < 0.20(0.16) GeV^/c^ in the one 
(three) prong search selects di-jet topologies. Due to the higher average mul- 
tiplicity, the 3h topology is more sensitive to variables related to the internal 
structure of Ty (Section 5.4) and H, thus requiring a tighter constraint. 

An additional check of the H structure is obtained by removing Ty and by 
computing the maximum transverse momentum QLep among all charged tracks 
within the hadronic system. As explained in Section 5.3.2, events originating 
from I'ei'^e) CC interactions are indeed characterized by large values of Qhep for 
the leading lepton. The condition Qf^^p < 5.0 GeV/c further suppresses topolo- 
gies where the unidentified leading electron (positron) is erroneously included 
in the hadronic system (Figure 4b). 



5.3 Lepton veto 



Due to the presence of a highly isolated track (leading lepton), unidentified 
CC interactions can fake hadronic r decays. As discussed above, the selection 
of a single leading charged track defines two different topologies in CC back- 
ground events, depending on whether the leading lepton is the chosen track 
(Figure 4a) or is included in the remaining hadronic system H (Figure 4b). 
In the latter case the transverse plane kinematics is significantly distorted by 
the selection of a random leading particle, thus reducing their effectiveness. 
Specific algorithms are then developed in order to tag the leading lepton. The 
selection can be divided in two steps: 
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a) 



b) 
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Fig. 4. Possible topologies for an unidentified ^^{ug) CC interaction after the se- 
lection of a leading charged track in the analysis: a) the primary lepton is chosen 
as the leading particle; b) the primary lepton is included in the hadronic system 
H. Similar considerations can be applied to the ry selection (l is replaced by ry )• 
However, for the I7, 27 and 3h samples the choice of the primary lepton as one of 
the Ty particles can produce additional topologies, as explained in Section 5.4 and 
shown in Figure 5. 

• First level veto (FLV). 

Lepton identification criteria are used to reject events containing clear pri- 
mary leptons of any charge. 

• Second level veto (SLV). 

A kinematic criterion based on event topology is used to tag a single lepton 
candidate, Iqq or Iqq , among all negatively charged primary tracks. Subse- 
quently, additional lepton identification criteria are imposed to the selected 
track. 

This method aims at avoiding tight lepton identification criteria. In partic- 
ular, it is also possible to analyze events outside the geometrical acceptance 
of the relevant subdetectors (Figure 1). The analyses described in Ref. [7] 
were instead based on tight identification conditions and rejected all events 
containing high pt tracks escaping the detector. This resulted in a significant 
efficiency loss. The kinematic tagging plays a crucial role in the veto algorithm 
because events in which the leading lepton is correctly tagged are efficiently 
rejected by kinematics (Section 5.5). 

5.3.1 Muon tagging and veto 

The FLV rejects events containing a positively identified by the re- 

construction algorithms or by the presence of residual hits and track segments 
associated to a charged track in both stations of the muon chambers. 

The muon tagging in the SLV is based on the likelihood function: 
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Table 1. Number of observed events and expected background as a function of cuts for all the hadronic DIS decay modes. The 
corresponding signal efficiency (%) is listed in the first column (uj- CC). The effect of the selection on the r+ control sample is also 
shown. The numbers are fully corrected with the data simulator and the last column shows the ratio between data and predicted 
background. The efficiencies refer to the inclusive h{mr^) and 3h{mr^) samples and do not include branching ratios. Overlap events are 
included in each topology. 



Sample 




V-r CC 


NC 






CC 


I'e 4- 


CC 


Tot Bkgnd 


Data 


Data/Bkgnd 


Charge 




- 


- 


+ 


- 


+ 


- 


+ 


- 


+ 


- 


+ 


- 


+ 


TV candidate 


O7 


43.7 


80739 


88663 


27632 


30852 


1612 


1426 


109883 


120941 


114012 


126013 


1.04 


1.03 




I7 


43.1 


96869 


103264 


26631 


29094 


1670 


1591 


123070 


133949 


126128 


137616 


1.02 


1.03 




27 


33.7 


79934 


80613 


21247 


26742 


1358 


1494 


102539 


107849 


105399 


109286 


1.03 


1.01 




3h 


44.4 


105524 


116237 


37874 


41406 


1612 


2452 


145060 


160095 


148420 


166094 


1.02 


1.03 


H structure 


O7 


36.3 


68383 


78517 


23987 


23420 


1305 


940 


93675 


102877 


96020 


105938 


1.02 


1.03 




I7 


36.5 


81682 


90164 


20347 


22476 


1233 


1002 


103262 


113642 


104609 


115531 


1.01 


1.02 




27 


27.7 


66257 


68570 


16533 


19187 


976 


917 


83766 


88674 


85102 


89317 


1.02 


1.01 




3h 


32.8 


88472 


98737 


28794 


31171 


1179 


1576 


118445 


131484 


116106 


134206 


0.98 


1.02 


Lepton veto 


O7 


23.7 


39319 


40077 


7422 


16299 


165 


389 


46906 


66766 


47673 


64986 


1.02 


0.97 




I7 


17.5 


32128 


32260 


3930 


5598 


164 


289 


36222 


38147 


36818 


39066 


1.02 


1.02 




27 


12.6 


23669 


22800 


2847 


4345 


120 


244 


26536 


27389 


27168 


27913 


1.02 


1.02 




3h 


13.2 


19701 


25712 


2196 


2566 


182 


381 


22079 


28669 


21929 


27889 


0.99 


0.97 


In A™ cut 


O7 






























I7 


9.5 


8668 


7888 


718 


1499 


69 


126 


9355 


9612 


9212 


9408 


0.98 


0.99 






7,1 


8702 


8048 


843 


1515 


59 


113 


9604 


9676 


9840 


9951 


1.02 


1,03 




■Ml 


4.!) 


24ii 


2:!66 


147 


219 


28 


45 


258(i 


26:UI 




2684 


i.lj:', 


1,02 


InA'^'^ cut 


O7 


9.7 


12020 


11464 


717 


684 


66 


88 


12792 


12236 


13471 


12101 


1.06 


0.99 




I7 


8.6 


7816 


7469 


565 


1127 


62 


107 


8443 


8693 


8346 


8696 


0.99 


1.00 




27 


6.2 


7864 


7609 


616 


1137 


61 


96 


8531 


8741 


8789 


9097 


1.03 


1.04 




3h 


4.2 


2079 


2026 


119 


173 


20 


43 


2219 


2241 


2263 


2243 


1.02 


1.00 




O7 


1.3 


6.7±2.1 


7.3±2.5 


0.4±0.4 


2.0±0.6 


0.7±0.1 


4.0±0.3 


7.8±2.2 


13.3±2.6 


10 


11 


1.28±0.54 


0.83±0.30 




I7 


1.4 


6.7±2.1 


8.9±3.1 


0.5±0.4 


3.6±0.9 


0.9±0.1 


4.6±0.4 


8.1±2.2 


17.0±3.2 


8 


17 


0.99±0.44 


1.00±0.31 




27 


0.8 


3.4±1.3 


7.9±1.7 


1.1±0.6 


3.4±0.9 


0.9±0.2 


3.3±0.4 


6.4±1.5 


14.6±2.0 


6 


16 


1.11±0.55 


1.10±0.31 




3h 


1.3 


2.8±1.4 


5.3±1.2 


1.5±0.6 


2.6±0.9 


0.6±0.2 


2.0±0.3 


4.9±1.6 


9.9±1.6 


3 


10 


0.61±0.40 


1.01±0.36 



tions with unidentified muons. These events are characterized, in general, by 
large values of y^j, and the effect of variables relating the hadronic system 
to the muon track is weaker. For this reason, Cy only includes variables di- 
rectly describing properties of the track being considered. The 2D correlations 
between variables used to construct Cy are shown in Figure 3 (right plots) . 

A likelihood ratio. In , is built from z/^ CC events surviving the FLV, as the 
ratio of £y for the true muon and for other tracks. In the three prong search 
Iqq is defined as the track with maximum InA^ among all the negatively 
charged primary tracks. In the one prong search the tagging efficiency can be 
further increased by considering the bias induced by the ry selection, which 
returns a single charged track. In this case the ry charged track is also tagged 
as Iqq if, for this track, InA^ > 0, otherwise /^c chosen as the track with 
maximum In A^ (the most likely muon) . 

The algorithm correctly identifies the true muon in 94%, 89%, 89% and 86% 
of CC surviving the FLV for the O7, I7, 27 and 3h topologies. In addition, 
the muons not identified by the algorithm have very low energy and are thus 
generally included in the hadronic system H (Figure 4b). This results in a 
NC-Iikc configuration which is efficiently rejected by kinematics against NC 
interactions (Section 5.5). 

Only events in which Iqq is compatible with being a minimum ionizing par- 
ticle in both ECAL and HCAL are rejected at this stage. The suppression 
of the remaining z/^ CC background is achieved by exploiting the complete 
event kinematics under the assumption that Iqq is the leading particle (see 
Appendix A), as explained in Section 5.5. 

5.3.2 Electron veto 

Although i/g (Pe ) CC events represent only a tiny fraction of the full neu- 
trino interactions (Section 3), z/g (Pg ) can be potentially more dangerous than 
z/^ CC events. This is due to the lower identification efficiency of electrons with 
respect to muons (limited angular acceptance of TRD, efficiency loss for tt re- 
jection, bremsstrahlung emission, etc.) and to the fact that this effect is not 
restricted to specific topologies. In addition, have a harder energy spectrum, 
because they originate mainly from K decays, and therefore the primary lep- 
ton is well isolated from the hadronic jet. However, in the hadronic r decays 
the presence of 7's from tt" (sometimes converting in the DC volume) can pro- 
duce electron-like signals. For this reason, the analysis does not use stringent 
constraints against fe(Pe) CC interactions, which would result in a significant 
loss of efficiency. 

The FLV is designed to reject two different classes of events. First, events con- 
taining a primary track which is positively identified as e~{e~^) by both TRD 
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a) 



b) 




|i(e) 




V 



V 



Fig. 5. Internal ry topologies (I7, 27 and 3h) for CC baekgrounds in which the 
leading lepton is part of the ry candidate: a) ^^{I'e) CC interaction with part of 
TV embedded in the hadronic jet; b) CC interaction with bremsstrahlung emission. 



and loose PRS requirements are rejected if, for this track, p^/p^ > 1.0 and 
> 1.0 GeV/c. These are mainly high pt e~[e^) in events with small missing 
transverse momentum. A second background source originates from z/e(z/e) CC 
interactions where no particle identification is available because the primary- 
electron (positron) emitted most of its energy as a hard bremsstrahlung 7 and 
missed the relevant subdetectors. Therefore, for each primary track (of either 
charge) not reaching the TRD, we search for potential bremsstrahlung secon- 
daries with direction tangential to the track. The presence of such additional 
7's or charged particles is then used as a kind of "electron tag" and the event 
is rejected if, for the fully reconstructed lepton (sum of primary track and 
bremsstrahlung secondaries), p^/Pt > ^-^ Pt > 1-0 GeV/c. 

Since, contrary to the muon case, the electron identification inefficiencies are 
not restricted to specific topologies, the SLV tagging criterion is simply based 
on In . A I^q candidate is then defined as the ry charged track in one-prong 
events and as the negative ry track with the larger momentum uncertainty 
in three-prong events. Events from CC interactions in which the primary 
electron is chosen as r daughter candidate are indeed dangerous because they 
are highly isolated and cannot be suppressed by studying the jet structure. 
This happens in 85%, 62%, 57% and 51% of all events for the O7, I7, 

27 and 3/1 topologies. These events are rejected if the Iqq track fulfills loose 
electron identification criteria based on TRD and on the combined PRS-ECAL 
information. 

As described in the following, the background in the most sensitive region of 
the analysis (Table 5) consists almost entirely of t'e(^e) CC interactions. This 
is partially explained by the kinematic approach used (Section 5.5). 
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Fig. 6. Left plots: correlations between kinematic variables used to construct In A^^ for MC events of type NC (a,c,c) and r Vrir 
decays with the selection of the right tt as leading particle (b,d,f). Right plots: correlations between kinematic variables for MC events 
of type CC with unidentified ji (a,c,e) and r Ut-tt decays with the chosen Iqq as leading particle (b,d,f). 



5-4 Structure of the ry candidate 



For the I7, 27 and 3h topologies the internal Ty structure provides further 
discrimination against backgrounds. The corresponding information, incorpo- 
rated into the £™ function, has been already used for the ry identification 
(Section 5.1), through the likelihood ratio between correct and random com- 
binations in signal events. 

Once the Ty selection is performed, background events can still have a ry 
internal structure inconsistent with r decay. In particular, the CC sample 
contains two specific topologies. The first one includes events where ry con- 
sists of the leading /i(e) and of additional particle(s) from the hadronic jet 
(Figure 5a). A second possibility arises when the leading electron from a Ue 
CC interaction undergoes hard bremsstrahlung radiation and ry consists of 
different pieces of the original electron (Figure 5b). Both cases can be sup- 
pressed by a constraint on the ry internal structure, which makes use of the 
presence of p or ai resonances in r decays. However, the ry internal structure 
is expected to be less effective against NC background because these events 
may contain genuine resonances inside the hadronic jet. 

A likelihood ratio, In A^^ , is then built as the ratio of the C^^ function between 
the true combination in signal events and the selected Ty in a weighted mixture 
of and CC interactions. A loose cut at InA^'^ > 0.5,0.0,0.0 is apphed 
for the I7, 27 and 3h topologies respectively. The higher charged multiplicity 
of the 3h topology partially correlates the ry internal structure to the jet 
structure (Section 5.2) for backgrounds and therefore a further binning is used 
along InA^^ in this case (Section 5.5). Since no kinematic constraint on any 
variable related to the ry structure has been previously imposed, the In A^'^ cut 
results in a significant event reduction for all samples (Table 1). However, this 
is due mainly to the rejection of events which do not contain Ty candidates 
compatible with being a genuine p or oi resonance. 



5.5 Global kinematics 

Background events can be divided into two categories with opposite kinematic 
configurations. In NC interactions the Ty candidate is embedded in the hadron 
jet (Figure 2a), and a large missing transverse momentum associated with the 
escaping neutrino is almost opposite to the direction of the hadronic system. 
On the other hand, the leading lepton in CC interactions is typically well- 
isolated and balances the momentum of the remaining hadronic system in the 
transverse plane (Figure 2c). The signal from hadronic t decays has interme- 
diate properties between these two extremes; the t decay neutrino introduces 
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Fig. 7. Scatter plot of In A*-^^ vs In A^^ for MC backgrounds from NC (left plot) and CC interactions passing the first level veto 
(middle plot) and simulated r — > i^rTr (O7) decays (right plot). The local bin densities are shown in logarithmic scale. The large box 
at the upper right corner, divided into four bins, defines the signal region, whereas the dotted lines show the control regions used to 
estimate the systematic uncertainties on CC (region A) and NC (region B) events. 



a modest amount of missing transverse momentum and the non-colinearity of 
T and Tycan reduce the isolation of Ty (Figure 2b). 

In order to optimize separately the rejection of each of the two opposite back- 
ground sources, we implement an event classification based on the use of two 
distinct likelihood functions. As discussed in the following, this procedure pro- 
vides the best description of the kinematic information and, at the same time, 
the possibility to distinguish further between NC and CC backgrounds. The 
first likehhood function is designed to separate the signal from NC interac- 
tions: 

C^^ = [[[(^pT,0„H],OTvhi:QT],PT^PT ] 

where the minimum opening angle 6^^ h. is sensitive to the internal structure 
of the hadronic system. Figure 6 (left plots) illustrates the discriminating 
power of some correlations between the kinematic variables used to construct 
jC^'^ . For each event, a hkelihood ratio InA'^*-' is computed as the ratio of 
the C^^ functions constructed from the true ry combination in signal events 
and from the selected ry combination in NC events respectively. Since C^'" 
is built by selecting ry as the leading particle, for the CC background a cut on 
suppresses the topology where the unidentified lepton has little energy 
and is part of the hadronic system H (NC-like configuration in Figure 7). As 
explained in Section 5.3.1, this reduces the mistagging probability of muons 
to a negligible level in the final signal region. 

The second function is designed to distinguish signal from CC events and is 
optimized, in particular, to reject CC interactions: 

= [[RQ^,Rp^,6^i],E^is,PT7MT] 

where the first part is similar to the function used for the muon tagging pro- 
cedure (Section 5.3.1). The function >C*-^*-^uses variables referring to the tagged 
muon Iqq and therefore, for events where the lepton is correctly tagged, 
it represents the actual CC kinematics (Figure 4). The potentially different 
choice of leading particle(s) provides, globally, four additional degrees of free- 
dom (including the Rq^ variable) as explained in Section 5.1, and thus further 
justifies the use of two distinct likelihood functions. The correlations between 
some kinematic variables included in are shown in Figure 6 (right plots) . 
For each event, a hkelihood ratio In A*-'*-' is computed as the ratio of the £^*-' 
functions constructed from the selected Iqq in signal events and from the true 
muon in MC CC interactions passing the FLV. 

Events are plotted in the plane [In A'^^ , In A*-'*-' ] , as shown in Figure 7 for sim- 
ulated signal (O7) and backgrounds. Two distinct populations, corresponding 
to the NC and CC backgrounds, are clearly visible, demonstrating that in- 
deed the two background sources are independent and have opposite overall 
kinematics. The signal region (the blind box), further subdivided into differ- 
ent bins, lies at large values of both hkelihood ratios. This corresponds to the 
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Fig. 8. Definition of the signal region: sensitivity as a function of In A^^ for the O7, I7, 27 and 3h (from left to right) topologies. The 
vertical lines correspond to the binning used in the analysis. The blind region starts at InA'^'^ > 5.5,6.0,5.5,4.0 for the O7, I7, 27 and 
3^ topologies respectively (left vertical lines in each plot). A cut InA*-"-^ > 3.5,0.0,0.0,0.0 is applied. 



selection of events where the ry candidate is isolated from the hadronic system 
and is not fully balanced in the transverse plane. 

For all topologies the analysis reaches its optimal sensitivity at low back- 
ground levels (Table 5). This is an indication that the kinematic rejection of 

background is optimal. The sensitivity curves of Figure 8 clearly show a fast 
rise after the minimum at large values of the likelihood ratios, associated to 
regions of decreasing signal efficiency and almost background-free. The pre- 
liminary quality cuts on the statistical significance of the main variables (Sec- 
tion 4) allows us to exploit such extreme values of the likelihood ratios. Due to 
the efficient kinematic suppression of the main background sources (NC and 
z/^CC), the very small residual background in these regions consists almost 
entirely of Vf. [pf. ) CC events. 

The choice of the binning for the signal region is performed before analyzing 
data events and is the one providing the best overall sensitivity to oscillations 
(Figure 8). This criterion allows an objective definition of bins. First, the 
highest bin is chosen such that its lower edge is as close as possible to the 
minimum of the curves of Figure 8, while consistent with a background of 
less than 0.5 events. Additional bins are then included in the signal region, 
up to the point where no further significant improvement of the combined 
sensitivity is obtained. The final binning follows the scheme of Figure 7 for O7 
and I7 events. In the 27 topology, due to the hmited statistics, only In A'^*-' is 
used to define the binning. For the 3h topology, the signal bins use also In A^'^ 
as a third dimension, since the higher multiplicity enhances the effect of the 
internal ry structure on background rejection. The last bin along In A^*^ is thus 
further divided into three bins along In A™ and In A*"^ respectively. 

Table 1 summarizes the analysis flow on the various samples. The signal region 
is defined as InA^^ > 3.5,0,0,0 and InA^^ > 5.5,6.0,5.5,4.0 for the O7, 
I7, 27 and 3h topologies respectively. The absence of a ry internal structure 
necessitates a tighter constraint on In A^*" for the O7 topology. In the 3h decays, 
due to the large ai mass, a smaller phase space is available for the final state 
Ur, thus increasing the effectiveness of the In A'^*-' cut. The final signal efficiency 
is similar for all topologies (Table 1). All signal bins are hsted in Table 2. 



6 Reliability of the background estimate 

Since a r signal would appear as a statistically significant excess of events 
inside the signal region, a crucial point for the analysis is the control of back- 
ground predictions in this region. This requires two distinct steps. First, cor- 
rections to the Monte Carlo are extracted from the data themselves. Then, 
data control samples are used to check final predictions and to evaluate the 
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Table 2 

Number of background and data events in the signal region. The corresponding 
N!^'^ and N^"^ , as defined in Sections 7.2 and 8.3, are listed in the last two columns. 
The 1/27 and O/I-27 topologies contain overlap events. The bins denoted by a star 
are considered as low background bins (Section 8.2). 



Analysis 




Bin# 


Tot bkgnd 


Data 








vrh{n^^) DIS 


O7 


I 


4.49 ± 1.50 


5 


454 


10.4 






O7 


II 


3.07 ± 1.17 


5 


345 


8.2 






O7 


III 


'-'•'-'^ -0.03 





288 


6.9 


* 




O7 


IV 


'-'•^^ -0.05 





1345 


31.1 


* 




I7 


I 


4.47 ± 1.58 


5 


283 


6.8 






I7 


II 


1.54 ±0.89 





244 


5.7 






I7 


III 


-0.04 





223 


5.7 


* 




I7 


IV 


-0.04 





1113 


26.6 


* 




27 


I 


2.57 ±0.91 


3 


318 


7.4 






27 


II 


0.66 ± 0.44 





175 


4.1 






27 


III 


0.49 ± 0.40 





82 


1.9 






27 


IV 


11+°-^° 

^■^^ -0.06 





211 


4.9 


* 




1/27 


I 


1.40 ±0.77 


2 


154 


3.7 






1/27 


II 


17 + 0™ 

-0.08 





124 


2.9 






1/27 


III 


20 + ° ™ 

•J-^U -0.06 


1 


707 


16.9 


* 




0/1-27 


IV 


14 + 0-70 

"•-■-^ -0.06 





1456 


34.2 


* 


z/^3/i(n7r°) DIS 


3/i 


I 


2.61 ±0.99 


2 


170 


4.0 






3/i 


II 


0.58 ±0.57 





139 


3.4 






Zh 


III 


0.86 ±0.57 





74 


1.7 






3h 


IV 


0.55 ± 0.59 


1 


309 


7.6 






3h 


V 


^■'^■^ -0.32 





675 


16.6 


* 



corresponding systematic uncertainties. 



6.1 Data simulator corrections 

The procedure to evaluate the backgrounds is based upon the data simulator 
method (Section 4.2), which imphes the apphcation of the full selection scheme 
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to appropriate data (DS) and simulated (MCS) samples. However, the over- 
all background is actually composed of three different physical sources: NC, 
^ni^fi) CC and fei^^e) CC interactions (Table 1). In order to obtain a precise 
background evaluation, these three categories are considered separately: 

• The NC correction is obtained from identified CC interactions, by the 
removal of the leading muon. This results in hadronic systems of different 

charge distribution than genuine NC events. This charge bias is corrected 
for by selecting ry of both positive and negative charge and by averaging 
the two results. 

• The i^nii^fj,) CC correction takes into account two distinct effects: the muon 
identification efficiency (veto) and the actual kinematic selection. The first 
correction factor is measured from a large sample (~ 6 x 10^) of muons 
originating from a nearby beam and crossing the NOMAD detector. The 
kinematic correction factor is evaluated from identified CC interactions in 
which the identification of the leading muon is ignored and, instead, the full 
T selection is performed. The topological bias related to unidentified muons 
(muon chamber acceptance, low muon momenta) has little effect since this 
class of events is efficiently rejected by the kinematic veto (Section 5.3.1) 
and is checked through an appropriate CC control sample (Section 6.2). 

• The z^e(^e) CC correction is obtained from events identified as CC in the 
T — > ePgi^r DIS search [7]. These events are then passed through the full 
event selection by ignoring the lepton identification. The electron identifi- 
cation (veto) used in the analysis, mainly based on TRD requirements, was 
independently checked using a sample of high-energy 5 rays produced in the 
DC volume by muons originating from the nearby beam. As explained in 
Section 5.3.2, unidentified Ve{ve) CC do not have specific topologies. There- 
fore, the described correction is adequate for all kinematic configurations. 

The total net correction factors, eDs/cMCSi to the number of background events 
computed from the Monte Carlo are 1.8 for O7, 2.0 for I7, 1.3 for 27 and 1.1 
for 3/i topologies respectively. These numbers are integrated over the whole 
signal region and are dominated by the data simulator corrections of the NC 
samples. In the most sensitive kinematic region, characterized by only a small 
residual CC background (Table 5), the corresponding correction factors are 
0.70 for O7, 0.91 for I7, 0.96 for 27 and 0.99 for 3/i topologies respectively. 

The data simulator corrections computed for Vr CC events are 0.86 for single 
prong and 0.99 for three prong events and are essentially independent of the 
values of the hkelihood ratios. 
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10 10 10 2 8 

In X'^^ 

Fig. 9. Cumulative InA^*^ distributions in control samples for the O7, I7, 27 and 3h topologies (from left to right): a) after cuts on 
the structure of the hadronic system H (Section 5.2); b) CC control samples complementary to the signal region; c) final distribution 
outside the signal region (difference between the previous two); d) wrong sign analysis. The histograms represent the individual NC and 

CC contributions to the total background. The vertical lines show the starting point of the signal region (solid) and of the last signal 
bin (dashed). The region to the right of the vertical solid lines in plots a) and b) is used to estimate the systematic uncertainties on the 
CC sample (region A of Figure 7) . 



6.2 Control samples 



The definition of appropriate control samples, which are needed to validate 
the final background predictions, must take into account two problems. First, 
the possibility to independently check each of the individual background con- 
tributions is desirable and requires, in turn, the ability to discriminate among 
them. In addition, the statistical uncertainty associated with each control sam- 
ple must be small in order to provide a significant test of systematic effects. 

The selection scheme based on two distinct likelihood functions provides by 
construction a separation between NC and CC backgrounds. Therefore, the 
required control samples can be defined in a natural way by selecting different 
regions in the plane of Figure 7. Since In A^*-' is built only from the CC 
sample, this latter background can be further partially decoupled from Vg 
CC interactions in the plane of Figure 7. This is important since in the low- 
background region essentially only fe^^e) CC events are present. However, for 
the sake of clarity, we will regroup backgrounds into NC and CC {v^ and Ve) 
interactions in the following, unless otherwise specified. 

The CC control sample is defined by analysing the projection along In A'^*" of 
events which fail at least one of the following selection criteria: 

• lepton veto; 

• In A™ cut; 

• InA*-^*^ cut. 

At large values of In A'^*-' this sample is completely dominated by CC interac- 
tions, as can be seen from Figure 9a-b. 

Conversely, the NC control sample is defined by analysing the projection along 
In A^*-' of the events which fail the In A'^*-' cut, after applying all the remaining 
cuts. Figure 10 shows the statistical significance of the individual NC and CC 
contributions to this sample (see also Figure 11). 

In addition, the distributions of InA'^*-' for the r"*" sample, where no signal 
is expected because of the small content of the beam (Section 3), are 
compared with data for each of the individual steps of Table 1. A similar check 
is performed in the r~ search for the initial selection and for events outside the 
signal region (Figure 9). Due to the charge bias of the lepton tagging procedure 
(Section 5.3) and, consequently, of kinematics (Section 5.5), the r"*" selection 
is less effective than the t~ search in rejecting backgrounds (Table 1). This 
gives the possibility to check background predictions with larger statistics. 

Data events are in good agreement with background predictions for all control 
samples, as shown in Table 1, Figure 9, Figure 10 and Figure 11. This gives 



28 




Fig. 10. Cumulative In A distributions of data (full circles) and background predictions (open squares) in NC control samples for 
the O7, I7, 27 and 3h (from left to right) topologies: a) without cuts on In A^*-' and In A^^ ; b) with a In A^^ cut reducing the overall 
background by more than a factor of 10 but without In A^'^ cut; c) with InA^'^ > 0.5,0.0,0.0 (signal selection) but without In A'^'-' cut. 
The histograms represent the individual NC and CC contributions to the total background. The vertical lines show the starting point 
of the signal region (solid) and of the last signal bin (dashed). The region to the right of the vertical solid lines is used to estimate the 
systematic uncertainties on the NC sample (region B of Figure 7). 
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Fig. 11. Cumulative In A distributions of data (full circles) and background predic- 
tions (open squares) for the I7, 27 and 3h (from left to right) topologies: a) without 
cut on InA*^*^ ; b) with In A'~'^ > 0.0, 0.0, 0.0 (signal selection). The histograms rep- 
resent the individual NC and CC contributions to the total background. The vertical 
lines show the starting point of the signal region (solid) and of the last signal bin 
for 3h events (dashed). The region to the right of the vertical solid lines is used to 
estimate the systematic uncertainties. 

confidence in the background estimation procedures and allows at the same 
time an evaluation of the systematic uncertainties. 



7 Systematic uncertainties 

The use of likelihood functions incorporating the full event topology provides 
a better estimate of systematic uncertainties with respect to a selection based 
on the application of a sequence of cuts. Moreover, the separation between 
CC and NC backgrounds gives additional high statistics control samples to 
constrain background predictions, as described in Section 6.2. This results, in 
turn, in a more precise background estimation. 

7. 1 Background 

The systematic uncertainties on the number of background events predicted 
inside the signal region can be divided into two contributions, related to the 
overall normalization and to the shape of the likelihood ratios. 
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The first term is estimated from a comparison of the total number of data and 
background events in both r+ and searches as a function of the selection 
criteria. This results in a r.m.s. of 1.9% from the differences in the last columns 
of Table 1. This value includes uncertainties related to fluxes, cross-sections 
and to the effect of the selection cuts for CC interactions, which define the 
normalization, and NC events, which are the dominant background sample. 
For the remaining components (different from ) there is an additional sys- 
tematic uncertainty coming from the integral fiux predictions [12]. An upper 
limit on this contributions is obtained from the average systematic uncertain- 
ties on the corresponding energy spectra, which amount to 5.1%, 8.5% and 
12.0% for i/e , P^and PgCC interactions respectively. 

The second term is evaluated from the two control regions shown in Figure 7 
for CC (region A) and NC (region B) backgrounds respectively. This is ob- 
tained by analyzing the corresponding control samples, defined in Section 6.2 
as the event projections along each of the two likelihood ratio axes (Figures 10 
and 9a-b). The systematic uncertainty is then estimated as the r.m.s. of the 
distribution of differences between data and predictions inside the signal re- 
gion of each control sample. This procedure is based on the assumption that 
in the region A (region B) a In A'^*-' (In A'^*-' ) cut does not introduce discrepan- 
cies between data and predictions which are dependent on the In A^^ (In A*-"" ) 
distribution. Due to the limited rejection factor of the In A^*" cut (Table 1), 
this condition is fulfilled by the CC control sample and the resulting number 
(5.5%) thus provides the final uncertainty on the CC background shape. 

The tighter constraint imposed on In A'"^*^ requires additional checks of the ef- 
fect of this variable on the In A*"*" distributions in region B. This is achieved by 
comparing the level of the observed agreement with different cuts on In A^*" , 
chosen in such a way as to reduce the overall background by more than a factor 
of 10 with respect to the initial value (Figure lOa-b). Since the Ty structure em- 
bedded in In A^'^ is almost independent of the global kinematics (apart from an 
overall boost factor in the laboratory frame), the cut on In A™ is not applied, 
in order to increase the available statistics and to have a consistent comparison 
of all topologies. The effect of the In A™ constraint is then checked separately, 
with no cut on In A'^'^ (Figure 10c). In addition, all the In A^'^ distributions for 
backgrounds are compared with data with and without the In A*"*" cut (Fig- 
ure 11). The final systematic uncertainty (r.m.s.) inside the signal region is 
estimated to be 4.1% from all these control samples. 

Table 3 summarizes the individual contributions to systematic uncertainties, 
which, added in quadrature, result in a total systematic uncertainty of 5.8% 
for z/^ CC, 10.3% for z/^ CC, 7.6% for z/e CC, 13.3% for CC, and 4.5% for NC 
respectively. However, the effect of this latter contribution is negligible since 
only the CC background is present in the most sensitive region of the analysis 
(Table 5). The overall net systematic uncertainty on the final background 
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Table 3 

Contributions to systematic uncertainties for background and signal events. 



Background 


Normalization 


1.9% 


Integral Ue/i^n 


5.1% 


Integral /v^ 


8.5% 


Integral Ve Iv^ 


12.0% 


CC InA shape in signal region 


5.5% 


NC In A shape in signal region 


4.1% 


Signal 


Normalization 


1.9% 




3.0% 


In A shape in signal region 


5.0% 



predictions for O7, I7, 27 and 3/i topologies is then 5.0%, 5.1%, 5.4% and 
5.5% at the starting point of the signal region and 10.0%, 9.2%, 9.2% and 
5.8% in the low background region. Given the large rejection factors of C(IO^) 
such a control of systematic uncertainties is noteworthy. 

1.2 Signal 

The probability of oscillation, Pqsc; is estimated as the ratio between the num- 
ber of observed r events and the maximal number of signal events expected if 
all incident f^had converted into Vr (Section 8.3). For the i/^ — > v^- oscillation, 
neglecting a term 0{P^^), this last quantity is defined by: 

where: 

• N'^^ is the observed number of z/^ CC interactions (Section 1). The number 
of z/^ CC interactions corresponding to the LM topologies are evaluated to 
be 11% of the total [7]. 

• Et and are the detection efficiencies for r signal events and i/^ CC events 
respectively, integrated over the incident spectrum. The cuts used to se- 
lect iV°^" and vary from channel to channel in order to reduce systematic 
uncertainties in the ratio / for that channel. 

• cTt/ C/i is the suppression factor of the Uj. cross section due to the difference 
between the r and masses, averaged over the incident i/^ spectrum. For 
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the spectrum used in this experiment [12] and for an energy-independent 
oscillation probability (corresponding to the large Am^ hypothesis), it is 
evaluated to be 0.48, 0.60 and 0.82 for the deep inelastic, resonance and 
quasi-elastic processes. The resulting average values for the DIS and LM 
analyses are 0.48 and 0.57 respectively. The event sample selected by the 
LM analyses contains significant fractions of DIS events satisfying the LM 
selection [7]. 

• Br is the branching ratio for the r decay channel under consideration. 

The systematic uncertainty on the overall normalization for the signal sample 
is common to the v^CG background (1.9%) and includes the effect of the 
individual selection criteria listed in Table 1. The maximum Data Simulator 
correction applied to the signal likelihood shapes, decoupled from the overall 
normalization, is 5.0%. We use this maximum correction as the systematic 
uncertainty on the signal shape. 

The uncertainty on the suppression factor Orl^n^ which will also contribute to 
the overall systematic uncertainty, is due to the uncertainties on the energy 
spectrum and on the structure functions used in the computation. It is esti- 
mated to be 3.0%. The uncertainty on the r branching ratios is negligible [25]. 

The final systematic uncertainty on the signal sample is 6.1% from the sum 
in quadrature of all the individual contributions summarized in Table 3. 



8 Results 

8.1 Analysis of the signal region 

After the choice of the selection criteria and the cross-check of the background 
predictions, we analyze data events falling in the signal region of the hadronic 
DIS channels (the blind box). Data events populate the various bins in a 
manner consistent with backgrounds, as summarized in Table 2. The overall 
integrals and shapes of the final In A distributions are in good agreement with 
background predictions (Figure 12). Therefore, no signal from oscillations is 
observed. 

8.2 Combined results 

The overall NOMAD results are obtained from the combination of the new 
analysis of the hadronic r decays with the Ur^i^e DIS and the LM analyses 
described in Ref. [7] (Appendix B). Table 4 fists all the individual contributions 
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Fig. 12. Histograms of In A for events passing all cuts except the cut on InA^ in a) O7, b) I7, c) 27 and d) 3/i, for data (points with 
statistical error bars), total backgrounds (shaded) and iV^'^ (Section 7.2), scaled down by a factor 20 (unshaded). The insets give, for 
each value of InA^*-' , the total number of events beyond that value, for data (dots) and expected background (squares); the encircled 
points are at the boundary of the signal regions. 



Table 4 

Summary of backgrounds and efficiencies for all the individual r searches. The 
columns labelled r~ summarize the observed number of candidate events (Obs.) 
and the corresponding predicted background (Tot Bkgnd) for the sum of all bins 
in the signal region. The columns labelled r"*" contain the equivalent numbers for 
the positive control sample. The corresponding r~ selection efficiencies (e,-), not 
including branching ratios, and the A^^"^ and N^'^ (Sections 7.2 and 8.3) are also 
listed. For the LM topologies the quoted e,- is the average efficiency for quasi- 
elastic and resonance events. The last column shows the overall sensitivity of each 
individual r search, based on the combination of all signal bins. 













r+ 










Analysis 


Obs. 


Tot Bkgnd 


Obs. 


Tot Bkgnd 


er (%) 






(xlO-4) 




DIS 


5 


r. O + 0.7 
-J-J -0.5 


9 


8.0 ± 2.4 


3.6 


4318 


88.0 


8.0 




DIS 


21 


19.5 ±3.5 


44 


44.9 ± 4.6 


2.2 


7522 


177.4 


4.0 




DIS 


3 


4.9 ± 1.5 


10 


9.9 ± 1.6 


1.3 


1367 


33.3 


22.2 




LM 


6 


5.4 ±0.9 


3 


2.2 ±0.5 


6.3 


864 


8.8 


55.2 




LM 


12 


11.9 ±2.9 


40 


44.1 ± 9.2 


1.9 


857 


16.7 


88.9 


^'T-3^(n7r°) 


LM 


5 


3.5 ± 1.2 


1 


2.2 ± 1.1 


2.0 


298 


5.2 


161.0 



Table 5 

Summary of background and data events in the low background bins. The corre- 
sponding N!^'^ and N^'^ , as defined in Sections 7.2 and 8.3, are listed in the last 
two columns. The quoted background consists mainly of CC events for the v-rei'e 
channel; of a mixture of and CC events for the Vrhimr^) channel and of z/^ 
CC events for the ^'T-3/i(n7r°) channel. The errors include contributions from the 
remaining sources. 



Analysis 




Bin # 


Tot bkgnd 


Data 




jyer 


UT-eUe DIS 




III 







680 


15.0 






VI 


0.16 ±0.08 





1481 


32.7 


(^vis < 12 GeV) 


II+III+VI 


0.27 ±0.13 





665 


8.7 


urhimr^) DIS 


O7 


III 


U-UO _o.03 





288 


6.9 




O7 


IV 


12 + 

^•-•-^ -0.05 





1345 


31.1 




I7 


III 


U-U' -0.04 





223 


5.7 




I7 


IV 


r, r,7 + 0.70 
U-U' -0.04 





1113 


26.6 




27 


IV 


11 +0-^0 

^•^^ -0.06 





211 


4.9 




1/27 


III 


^•^^ -0.06 


1 


707 


16.9 




0/1-27 


IV 


•J-^^ -0.06 





1456 


34.2 


i^r2h{mr^) DIS 


3h 


V 


"J-'J^ -0.32 





675 


16.6 


Total 






J^-"^ -0.39 


1 


8844 


199.3 
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Fig. 13. Data event (run 14162, event 25050) falling in the low background signal 
region of Table 5. The event is classified as r ^ pv^- (bin I/27 III) decay candidate: 
InA^'^ = 2.1, InA^^ = 9.4 and InA^*"" = 7.5. Solid lines represent the reconstructed 
charged tracks (open triangles show their extrapolation) whereas dashed lines rep- 
resent the ECAL neutral clusters. The energy deposition of individual ECAL cells 
is also shown by the shaded bins on the right. 

and the corresponding sensitivities. The unified analysis of Urh{mT^) topologies 
described in the present paper is the most sensitive appearance search in 
NOMAD. 

As discussed in Section 5.5, the strength of this search lies in the possibility 
to define, through kinematic constraints, large regions characterized by the 
expectation of low background. A precise control of background predictions 
using the data themselves is also crucial. The major low background (< 0.5 
events) bins are summarized in Table 5. Overall, these bins contribute about 
75% of the total sensitivity of the experiment. A single event is observed within 
this kinematic region. This event, shown in Figure 13, is classified as a r — > v^-p 
candidate and has large values of all likelihood ratios (Figure 12b-c). 

The overall systematic uncertainties for the hadronic DIS channels are given 
in Section 7. For the remaining modes the estimated systematic uncertainties 
are 20% and 10% on backgrounds and A^^"^ respectively [7]. 

8. 3 Evaluation of confidence regions 

The final result of the measurement is expressed as a frequentist confidence 
interval [27] which accounts for the fact that each r decay mode and signal 
bin (Tables 2 and B.l) may have a different signal to background ratio. The 
acceptance region of Ref. [27] therefore becomes mult i- dimensional to contain 
each of the separate measurements. The procedure follows the prescription of 
Ref. [25]. This computation [26] takes into account the number of observed 
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signal events, the expected background and its uncertainty, and the maximal 
number of the expected signal events. 

The resulting 90% C.L. upper limit on the — > u-,- oscillation probability is: 

Posc{^^^. ^^r)< 1-63 x 10"^ (2) 

Under a two-neutrino family formalism this corresponds to sin^ 2^^^ < 3.3 x 
10~^ for large Am^ and to the exclusion region in the Am^ — sin^ 26 plane 
shown in Figure 14. The result is significantly more stringent than the pre- 
viously published limits [7,20-23]. The sensitivity [27] of the experiment is 
Pose = 2.5 X 10""^; this is higher than the quoted confidence limit, since the 
number of observed events is smaller than the estimated background. In the 
absence of signal events, the probability to obtain an upper limit of 1.63 x 10~^ 
or lower is 37% (Figure 15). This result matches the design sensitivity of the 
experiment (Pqsc = 1-9 x 10~^). 

In the context of a two-flavour approximation, we can reinterpret the result 
in terms of — > ly-r oscillations, by assuming that any observed u^. signal 
should be due to oscillations from the small z/g component of the beam [28]. 
The corresponding maximal number of signal events, N!f'^ , is then obtained 
from N^'^ by reweighting the signal events using the i^e to flux ratio. This 
procedure introduces a further systematic uncertainty of 4.7% on the values 
of N^'^ , related to flux predictions [12]. The resulting 90% C.L. upper limit 
on the — > Vt oscillation probability is then: 

PoU^e ^Vr) < 0.74 X 10-2 (3) 

corresponding to sin^ 29^^- < 1.5 x 10^^ for large Am^. The exclusion region in 
the Am^ — sin^ 29 plane is also shown in Figure 14. The z/g — > Ur sensitivity 
is Pose = 1.1 X 10^2 and the probability to obtain an upper limit of 0.74 x 10"^ 
or lower is 39% (Figure 15). Both the sensitivity and the probability to obtain 
the result (goodness-of-flt) are an essential part of the limits themselves, as 
pointed out in Ref. [25] [27]. 

The results from the r appearance search also exclude effective couplings of 
1/^ or Ug with the r lepton, which are equivalent to the oscillation probabilities 
at large Am^. In particular, this information is required in order to relate 
the recent observation of r production by the DONUT collaboration [29] to 
the presence of z/^ in the beam. At 99% C.L. the NOMAD data limit such 
couplings to P^r < 4.4 x 10"'^ and Per < 2.0 x 10"^. 
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Fig. 14. Contours outlining a 90 % CL region in the Am^ — sin^ 26 plane for the two-family oscillation scenario. The NOMAD 
(left) and Vt (right) curves are shown as solid lines, together with the limits published by other experiments [20-24] 
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Fig. 15. Histograms of the upper limits obtained, in the absence of signal events, 
for 500 simulated experiments with the same NOMAD expected background [27]. 

The averages correspond to the quoted sensitivities while the arrows arc the actual 
upper limits obtained from the data in the Ufj, — > Vr (left) and — > i't (right) 
searches. 



9 Conclusions 



The analysis of the full NOMAD data sample gives no evidence for z/^ ap- 
pearance. In the two-family oscillation formalism this result excludes a region 
of the — > Vr oscillation parameters which limits sin^ 20^,- at high Am^ to 
values smaller than 3.3 x 10~^ at 90% C.L., and Am^ to values smaller than 
Am^ < 0.7 eV^/c^ at sin^ '29^r = 1- The corresponding excluded region at 90% 
C.L. for the — > Ut oscillation parameters includes sin^2^eT < 1-5 x 10~^ at 
large Am^ and Am^ < 5.9 eV^/c^ at sin^ 2^eT = 1- Our sensitivity to oscilla- 
tions is not limited by backgrounds, but is essentially defined by the available 
statistics. 



The NOMAD experiment has explored neutrino oscillations down to proba- 
bilities which are more than one order of magnitude smaller than limits set by 

the previous generation of experiments. For the first time, a purely kinematic 
approach has been applied to the detection of CC interactions. Our final 
results demonstrate that this approach has developed into a mature technique, 
providing a precise control of backgrounds from the data themselves. 
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Structure of H system 



Structure of Ty candidate 
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Global kinematics 
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Fig. A.l. Definition of the NOMAD kinematics for a i>r CC event. 
A Kinematic Vciriables 

The event kinematics is based on a set of global variables which describe 
the general properties of the two momentum vectors of a leading particle and 
of the hadronic system H recoiling against it, in the laboratory frame. This 
convention is closely related to the kinematics of CC interactions. However, 
in general, any track (or system of tracks) can be chosen as leading particle. 

Invariance with respect to an arbitrary rotation in the plane transverse to 
the beam direction means that an event can in fact be fully described by five 
degrees of freedom (Figure A.l): three in the transverse plane {x,y) and two 
along the beam direction (z) . 

In the selection of Uj. CC interactions, the leading particle consists of the vi- 
sible T decay product (s) Ty . The following kinematic variables can then be 
computed (Figure A.l): 

• Pvis, the total visible momentum of the event. This is computed by summing 
the momenta of all primary charged particles, neutral secondary vertices and 
neutral ECAL clusters. 

• Eyis , the total visible energy of the event. 

• and p^, the total momentum of the visible tau decay product (s) and 
of the associated hadronic system respectively, such that p"^^ + p^ = Pvis- 

• ysj, the ratio between and the total visible energy. 

• p^ and p-^, the components of p''^ and p^ perpendicular to the neutrino 
beam direction. 
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• Pj^, defined as — (p^^^ + Pt) ^^^^ interpreted as a measurement of the 
"missing" transverse momentum due to the neutrino(s) from r decay. 

• Mt, the transverse mass, given by = 4p^p^ sin^( 07^^/2) where (prym 
is the angle between and p^, when assuming massless decay product (s). 
For T events, Mt < r mass, up to detector resolution and Fermi motion 
effects. 

• Rpj., the ratio of the transverse momentum p!^ and the missing transverse 
momentum p^, 

• Qt, the component of p'^^ perpendicular to the total visible momentum 
vector (including Ty ) . 

• Qlcp , the component of a charged ])article momentum perpendicular to the 
total momentum of the rest of the event. 

• 9^H, the angle between the neutrino beam direction and the hadronic sys- 
tem. 

• 9^.ry , the angle between the neutrino beam direction and the Ty momentum 

vector. 

• 6uT, the angle between the neutrino beam direction and the total visible 
momentum vector of the event. 

• (^ryH, the angle between the hadronic system and Ty . 

For the lepton tagging algorithms (Section 5.3) and the kinematic rejection of 
CC interactions (Section 5.5) the lepton candidate track Iq^ is chosen as lead- 
ing particle in the computation of all the above variables (e.g. 9i,ry becomes 
Oi,i etc.). 

The following variables, partially related to the previous ones, incorporate 
information from the internal structure of the hadronic system H: 

• the minimum angle between Ty and any other primary track hi in 
the event. 

• {Qt)h, average computed among all the charged tracks of the remaining 
hadronic system H, after the exclusion of a particular track. This variable 
measures the transverse size of the hadronic system. 

• {Qt)t, average Q"^ computed among all the charged tracks in the event. 

• Rqj., ratio between the transverse size of the hadronic system, {Qt)h, and 
that of the full event, {Qt)t- This variable is sensitive to the isolation of 
the particle(s) not included in the hadronic system H. 

• Arryhi = \J {^V-TMhiY + {^4>Tyhi)'^, the minimum invariant opening cone 
between Ty and any other primary track hi in the event. This combines the 
differences of the corresponding angles in the transverse plane and of the 
pseudo-rapidity r] — — lntan(^/2). 

In addition, variables describing the internal structure of the candidate Ty are 
used, where applicable, in order to increase background rejection: 
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• Mp, invariant mass of a tt^tt*^ combination. 

• MpO, invariant mass of a n^n^ combination. 

• MttO , invariant mass of a 77 combination. 

• Mai ) invariant mass of a 7r~7r''"7r~ combination. 

• ^tt-ttO; opening angle between a 7r~ and a 7r°. The 7r° momentum can be 
obtained, in turn, from a single ECAL cluster (I7) or from the sum of two 
separate ECAL clusters (27). 

• 9^^, opening angle between two 7's. 

• ^TT+TT-; opening angle between a tt"*" and a 7r~. 

• ^^-TT", opening angle between two distinct 7r~'s. 

• E^o, energy of a tt^ obtained from a single ECAL cluster (I7) or from the 
sum of two separate ECAL clusters (27). 

• E^^^, maximum energy between two different 7's used to reconstruct a n^. 

• Etj.+ , energy of a tt"*". 

• EpO, energy of a tt'^tt" combination. 

where the track charges refer to the r~ selection (opposite for the r"*" selection) . 



B Additional topologies 

In addition to the analysis of the hadronic DIS channels described in this pa- 
per, the Ut appearance search in NOMAD includes the analysis of the — > 
e~Pe DIS decays and of the LM topologies from Ref. [7]. The corresponding 
values of A^^"^ and N!^'^ have been updated according to the most recent beam 
predictions [12]. Table B.l summarizes all the relevant numbers for these sam- 
ples, which supersede the ones quoted in Ref. [7]. 

In principle, NOMAD is sensitive to all leptonic and hadronic r decay chan- 
nels. However, the r~ fi'u^UT- decay channel is dominated by the z/^ CC 
background where the primary muon is positively identified in the detector. 
In view of the intrinsic difficulty of the evaluation of the systematic uncer- 
tainty on this background, due to the absence of a suitable control sample, 
this channel is not directly used for the Uj- search. 
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Table B.l 

Number of background and data events in the signal region for the t~ 
e~v^ Ut DIS and LM topologies [7]. The corresponding Nif''^ and N^"^ , as defined in 
Sections 7.2 and 8.3, are listed in the last two columns. The bins denoted by a star 
are considered as low background bins (Section 8.2). 



Analysis 


Bin# 


Tot bkgnd 


Data 






T~ — > p~D^ ij^ F)TS 


I 


85 + °?2 


2 


143 


2.8 


(^vis > 12 GeV) 


II 


^•^0 _ 0.12 


1 


136 


2.9 




III 


n 1 a +0.18 





680 


15.0 * 




IV 


1.85 ±0.22 


2 


554 


14.0 




V 


0.78 ±0.15 





406 


9.1 




VI 


0.16 ± 0.08 





1481 


32.7 * 


(^^vis < 12 GeV) 


I+IV+V 


0.77 ± 0.26 





253 


2.8 




n+ni+vi 


0.27 ± 0.1.3 





GG5 


8.7 * 




I 


3.09 ± 0.67 


3 


282 


2.9 




II 


1.50 ±0.41 


2 


286 


2.9 




III 


0.82 ±0.41 


1 


296 


3.0 


T ^ h{mT^) LM 


p 


5.2 ± 1.8 


7 


480 


8.9 


LM 


h 


6.7 ±2.3 


5 


377 


7.8 


T ^ 3/i(n7r°) LM 




3.5 ± 1.2 


5 


298 


5.2 
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